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We study an extended periodic Anderson model with the Coulomb interaction U c f between conduction and / elec- 
trons by the Gutzwiller method. The crossovers between the Kondo, intermediate-valence, and almost empty /-electron 
regimes become sharper with U c /, and for a sufficiently large U c f, become first-order phase transitions. In the Kondo 
regime, a large enhancement in the effective mass occurs as in the ordinary periodic Anderson model without U c f- In 
addition, we find that a large mass enhancement also occurs in the intermediate- valence regime by the effect of U c f. 

KEYWORDS: mass enhancement, Gutzwiller approximation, extended periodic Anderson model, valence transi- 
tion, valence fluctuations, heavy-fermion superconductivity 



1. Introduction 

In rare-earth and actinide compounds, several interesting 
phenomena, such as magnetism, heavy-fermion phenomena, 
and superconductivity, occur owing to the interplay of the 
strong Coulomb interaction U between / electrons and the 
hybridization V between the localized /-orbital and conduc- 
tion band. 

Among such phenomena, heavy-fermion superconductivity 
has been one of the central issues in /-electron physics after 
the discovery of the superconductivity in CeCu2Si2.^ In the 
heavy-fermion systems, the conventional, phonon-mediated, 
s-wave superconductivity is hardly realized owing to the 
strong onsite Coulomb interaction. Then, pairing mechanisms 
other than the phonon-mediated mechanism have been dis- 
cussed. The magnetic-fluctuation-mediated superconducting 
mechanism may be common in heavy-fermion superconduc- 
tors, since superconductivity is realized near a magnetic quan- 
tum critical point in many compounds. 

However, some heavy-fermion superconductors are diffi- 
cult to understand solely by the magnetic fluctuation scenario. 
For example, the superconducting transition temperatures un- 
der pressure in CeGi2Si2 2 ' and CeGi2Ge2, 3 ' become maxi- 
mum far away from the magnetic quantum critical points. In 
CeCu2Sii.8Geo.2, 4) the superconducting region splits into two 
regions: the low-pressure region close to the magnetic critical 
point and the high-pressure region away from the magnetic 
critical point. To explain the high-pressure superconducting 
phase, a valence fluctuation scenario is proposed 5-8 ' and the 
importance of the Coulomb interaction U c f between conduc- 
tion and / electrons has been discussed in addition to U and 
V. 

Valence fluctuations are expected from the rapid change in 
the valence of an / ion in these compounds under pressure, 
which is suggested from the behavior of the effective mass m*. 
In the periodic Anderson model (PAM), m* and the number of 
electrons, n f, in the / orbital per site follow the relation 9 ' 10) 

m* 1-W//2 

— = - — — , (1) 

m 1 — rif 

where m is the free-electron mass. This relation is derived 
using the Gutzwiller method for U — > oo. Thus, a large 
change in m* indicates a large change in tif. For CeCu2Si2 



and CeCu2Ge2, m* is experimentally deduced from specific 
heat measurements or the temperature dependence of elec- 
trical resistivity, and it is found that m* decreases rapidly at 
approximately the pressure where the superconducting transi- 
tion temperature becomes maximum. 11 ' 12 ' Thus, these obser- 
vations indicate that a sharp valence change or large valence 
fluctuations play important roles in superconductivity. 

However, there are problems with the relation between the 
effective mass and valence. First, eq. (1) is derived for an or- 
dinary PAM, which does not show a sharp valence change. 
Thus, we cannot naively apply eq. (1) to a system with large 
valence fluctuations. Second, the effective mass has a peak 
in CeCu2Si2 under pressure before the superconducting tran- 
sition temperature becomes maximum, that is, the effective 
mass varies nonmonotonically under pressure. 12 ' This non- 
monotonic variation in the effective mass cannot be explained 
by eq. (1), since the n/ of a Ce ion is expected to decrease 
monotonically under pressure for the following reasons. A 
positively charged Ce ion should be surrounded by negative 
charges. When these negative charges get close to a Ce ion un- 
der pressure, the / level ej of Ce is lifted. Under pressure, the 
overlap between the wave functions of the / orbital and the 
conduction band increases, and V increases. Both the effects 
of pressure on e/ and V result in a decrease in n f. Indeed, a 
monotonic decrease in n j underpressure has been observed in 
CeCu2Si2 by an X-ray absorption experiment recently. 13 ' We 
also note that, in CeCu2Ge2, the pressure dependence of the 
effective mass has a shoulder structure before the supercon- 
ducting transition temperature becomes maximum. 11 ' This 
shoulder structure may also become a peak as in CeCu2Si2, if 
we subtract the contributions of magnetic fluctuations, which 
are large in the low-pressure region. 

The peak structures in the effective mass may be explained 
by considering a combined effect of the renormalization de- 
scribed by eq. (1) and valence fluctuations. 12 ' However, the 
applicability of eq. (1) to a system with large valence fluctu- 
ations is not justified. Thus, it is an interesting problem how 
eq. (1) can be extended to a model that shows a sharp valence 
change. A study of such a problem will be helpful to under- 
stand superconductivity in CeCu2Si2 and CeCu2Ge2 by the 
valence fluctuation scenario. 

Note that heavy-fermion behaviors in Q'-YbAlB4 and /?- 
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YbAlB4 14) are also difficult to explain by eq. (1). The va- 
lences of the Yb ion are +2.73 for a-YbAlB4 and +2.75 for/3- 
YbAlB4, 15) that is, the hole numbers of Yb ions are n/ = 0.73 
and 0.75 for a-YbAlB4 and y8-YbAIB4, respectively. For these 
values of rif, i.e., much less than unity, we cannot expect 
heavy-fermion behavior from eq. (1). 

In this research, we study an extended periodic Ander- 
son model (EPAM) with the Coulomb interaction U c j be- 
tween conduction and / electrons by the Gutzwiller method. 
We extend the Gutzwiller method for the PAM developed by 
Fazekas and Brandow 10) to the present model. Then, we in- 
vestigate the effect of U c f on the effective mass. Although the 
EPAM has been studied as a typical model for valence transi- 
tion 1 ^ and has been investigated by some modern techniques 
in recent years 6-8 ' 17 ~ 19) after the proposal of the valence fluc- 
tuation scenario for superconductivity, the effect of U c j on the 
mass enhancement has not been clarified well. Some of the 
results have already been reported in our previous paper; 20) 
here, we report the details of the method and also add new 
results. 

2. Formulation 

The EPAM is given by 16) 

H =Tj £kc k'T ck < T+e fTj n f r ' r 



ka 



kcr 



(2) 



tlcratlfrcr' » 



where c k(r and faa are the annihilation operators of conduc- 
tion and / electrons, respectively, with the momentum k and 
the spin cr. n crcr and n j riT are the number operators at site r 
with <x of the conduction and / electrons, respectively, is 
the kinetic energy of the conduction electron. We have not 
taken orbital degrees of freedom into consideration. This sim- 
plification may be justified for a system with tetragonal sym- 
metry such as CeCu2Si2 and CeCu2Ge2 and for a system with 
orthorhombic symmetry such as a-YbAlB4 and y6-YbAlB4, 
since the crystalline electric field ground states of / elec- 
trons are Kramers doublets in these systems. In the following, 
we set the energy level of the conduction band as the origin 
of energy, i.e., Y^k e k = 0. We set U — > oo, since the on- 
site Coulomb interaction between well-localized / electrons 
is large. 

We consider the variational wave function given by 



where 



W) = PffPcM), 



P ff = fl [1 _n M"/<-|] 



(3) 



(4) 



excludes the double occupancy of / electrons at the same site, 
and 



Pcf = f~[[l - (1 - g)n cr o-nf rC r>] 



(5) 



rcrtr' 



is introduced to deal with the onsite correlation between con- 
duction and / electrons. 17) g is a variational parameter. The 



one-electron part of the wave function is given by 



(6) 



k<kj:,(T 



where kp is the Fermi momentum for the free conduction 
band without / electrons, |0) denotes vacuum, and a(k) is de- 
termined variationally. Here, we have assumed that the total 
number n of electrons per site is less than 2. 

In the method by Fazekas and Brandow, 10) the /-electron 
state is expanded in the basis state in real space, since it is 
convenient to deal with the projection operator Pff. In the 
present study, we consider P c f in addition to Pff, and thus 
we also expand the conduction-electron state in real space. 
This is the main difference from the method of Fazekas and 
Brandow. The creation operator is expanded as 



bl 



r bl = Yj<Pk{r)bl, 



(7) 



where L is the number of lattice sites and b denotes c a or fa-. 
Then, the basis state in momentum space is expanded as 

n„ 

nKi°> 



i{* (&) }> 



N b 



{rC>) 

The determinant is defined as 



= 2det[^,(r^)]f]fct.| 0> 

{!<»} i=\ 



(8) 



det[y>jt(r)] 



fkiin) 
<Pk 2 (ri) 



(9) 



The basis state including both c and / electrons is given by 

N ar N f<r 

(10) 



cr y=l ' i=\ ' 

where N ccr and Nf a are the numbers of conduction and / elec- 
trons, respectively, with spin cr. The total number of spin- 
cr electrons, No- = N ccr + TV/o-, should be fixed. Here, we 
have introduced the notation {r (e) } = {r (eT) , r (ci) ) and {r v) } = 
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r (/T) 5 Then, the one-particle part is expanded as and 

\4>)=Y\ J] ( _l)Perm((^}) 

f Ar /<r ,{tV")} 

7=1 '' 1=1 

* W/,.^") (11) 



xf]a(^) 

i=l 

x J] det[^,(r (OT) )]det[^ (/ „(r (/cr) )] 
x |{r (e) }{r (/) }>, 

where (-l)P erm (l* (/ ')) is the sign due to the fermion anti- 
commutation relation. In the summation, we should keep 
k (c<r) >jfc (Ar) < fep and {A . (OT)) n {k(fcr)} = % Then the 

tion in eq. (3) is carried out by restricting the summation in 
eq. (1 1) with the condition {r^Jnfr^} = and by multiply- 
ing each term by g D . D = D c f + D c j, where D C(T is the number 
in the set {r (c<T) }n{r ( ^}. That is, D is the number of interacting 
electron pairs through U c f. In the following formulation, we 
impose these restrictions without mentioning them explicitly. 

Then, we evaluate the normalization factor (tfr\tfr) using the 
approximation introduced in Appendix A. By using eq. (A-4), 
we obtain 

_ N fr Nfi 

z zn^on^) 

N fr N n {k W } i=\ 7=1 

x Yj g 2D |det[^ (£ ,,(r (cT) )]| 2 |det[ W /T,(r (/T) )]| 2 
M 

xldett^^Ml^detf^i)^)]! 2 , 

where = {Jk^.t^} and {r} = {r (c \r ( J">}. By applying 
eq. (A-3), we further approximate (i/f\i/f) and obtain 

(m* Z fxwfX'NfO 

N f iN n D cr D cl 

xY(Nf,N er ,D er )Y{N f ,N cl ,D cl ) (13) 
xZWfiWNfi), 
where tV/ = tV/t + tV/j, 



(12) 



X(N fr ,N fl ) = 



(14) 



7V f Cn C/y _£> 

YiNf^Dar) = / (15) 



jfcC/^j i=l 

= Z-p|-Zf- ln ^H 

(. i=l J 



(16) 



= exp[-F(/V /cr )]. 

Z(Nfo-) is the partition function of the canonical ensemble for 
the system with the dispersion = - In a 2 (£) with the con- 
straint k < k F at temperature 1 . F(Nfa-) is the free energy of 
this fictitious system. Then, by using the Stirling formula, we 
rewrite the normalization factor as 



<</#> - eX Pi L f( n fh n fU d ch d cl)l 

NnN fl D c ,D ci 



(17) 



where n/o- = Nfa-/L and <f co - = D C(T /L. In the summation, the 
most important terms should satisfy 

df(n fh n n ,d ch d ci ) 



dd a 



= 0, 



and 



df(n fh n fl ,d ch d cl ) 



dn fcr 

From eq. (18), we obtain 

n d aT (\-n f -n c 



= 0. 



8 



+ d ca ) 



(18) 



(19) 



(20) 



(n f - d CIT )(n co - - d c(T ) 

where n ccr = n lT -nf a with n a = No-/L. This is the same form 
as that in the Hubbard model, 21) if we regard n c<T as n^, rif 
as and d aT as d H , where n^i and c/ H are the numbers of 
cr-spin electrons and doubly occupied sites per lattice site, re- 
spectively, in the Hubbard model, and cr denotes the opposite 
spin of cr. From eq. (19), we obtain 



n 2 f (n co - ■ 



■ d ca ){\ - n ca ){\ - rif - n c£r + d cd ) 



(21) 



(1 - n ftT ){\ -n f )n CC r(n f - d c<T )(n f - d c& ) ' 

where fi is the chemical potential for the fictitious system de 
fined as 

dF(N) 



dN 



N=N fa 



(22) 



In the following, we assume a paramagnetic state, i.e., rif a = 
rif/2, n CCT = n c /2 = (n - «/)/2, and d ccr = d/2, and optimize 
the wave function so that it has the lowest energy. In the fol- 
lowing, we regard d as a variational parameter instead of g by 
using eq. (20). From the definition of the chemical potential in 
the grand canonical ensemble, the following equation should 
be satisfied 



B/ /2 = iy — -4 

L U 1 + e ' lna ' 



(k)-fi(n f /2) 



= l -y 

where we have introduced q 



a 2 (k) 



(23) 



■a 2 {ky 



If we set g = 1, that is, if we ignore the correlation be- 
tween the conduction and / electrons, we obtain d = n c rij 
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and q~ l = (1 - «//2)/(l - nf), which is the renormalization 
factor given in eq. (1). Our theory is reduced to the previous 
Gutzwiller method for the PAM by setting g = 1. This can 
also be checked for other quantities, such as renormalization 
factors, which will be derived in the following. 

Next, we evaluate the kinetic energy. The effect of the an- 
nihilation operator on the variational wave function is written 
as 

N tf N fl {r){ftW) 



xY]a(k^)Y]a(k ( f l) ) 



(24) 



x det (r,) [^(c T )(r (cT) )] det[^«/T,(r (/T) )] 
x det[^ W) (r (ei) )] det[<p km (r m) )] 



>T) ifp(c|)if_(/)i 



where 



det (r,) [^<c T ,(r (cT) )] 
^,, T »(r') ^(r^) 



(25) 



Thus, we need to introduce another approximation to evaluate 
the determinant eq. (25). By using eq. (A- 10), we obtain, for 
r + r', 

N fr N n D cr D cl 



xY(N f ,N ch D ct )Y(N f ,N cl ,D cl ) 



/=1 ' 

The renormalization factor q c(T is given by 
1 



(26) 



n c ,r(l -n ca ) 
X 



-^(n co - - d ccr )(l -rif -n aT + d ca ) 
+ ^d ca -(nf - d co -) . 



(27) 



q ca - has the same form as the renormalization factor q^ in the 
Hubbard model 21 ' as in the case of the Gutzwiller parameter 
g. From eq. (26), we can evaluate the momentum distribu- 
tion function n c (k) = n c] {k) = <4 T C *T> = {^\c[^c k] \tj/)l{^). 

In this evaluation, we need to calculate Z{*OT} Il^i a 2 (k^) 
with the restriction k t {k ( ^}. We can accomplish it with the 
aid of eq. (B-4). The result is 

/ j \ ((1 -?ctKt + An c( k ) for k<k F 
n C ](k) = \ , (28) 

(l-<3Vr)»cT forfe>fe F 



where 



An c (k) = q c 



a 2 (k) 



(29) 



In a similar way, we obtain rif(k) = n/i(k) = {f k ^fk-\) as 
1(1 - qfi)nfi + Arif(k) f 'or k < k-p 



n f i(k) = 



(1 -q f t)n n 



for k > kp 



(30) 



where 



a 2 (k) 



The renormalization factor for an / electron is 

r, - r, - 1 _ U f JfiT) Jcl) 

If ~ If* ~ If If ' 

1 flfcr 



(31) 



(32) 



where 



n f (l-n f ) 

X [ ^j(n f - d ccr ){\ -n f -n CIT + d CIT ) 
+ ^dcAn c(T - d co -)j ■ 



(33) 



q^ has the same form as q\\ in the Hubbard model, 21 ' if we 
regard «/ as n", n aT as n", and d ccr as d H . We can also evaluate 
the mixing term 



/ t * \ 1 V a( ^> 



(34) 



where the renormalization factor is given by 

(n f - d ca -)(n f - d ca -) 



q c f - qcfo- 



n 2 (l - n aT ) 



1 + 



(n f - d C o-)(l -n f - Ilea + d C cr) 



(35) 



Then, the expectation value of energy e - (H)/L per site is 
given by 

2 



k<kp 

2 
L 



2 ri fe- h)a 2 (k) - 2V x a(k) 



(36) 



z 



4- 1 + a 2 (£) 



where = q c e k and Vi = q c /V. We minimize the expectation 
value of energy with respect to the variational parameters a(k) 
and d. From de/da(k) = 0, we obtain 



a(k) = 



2Vi 



e f - h + 



P f -e k ) 2 +AV 2 2 



(37) 
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where V2 = yfqVi- The renormalized /-level ey should satisfy (a) n?- regime 



dq 
dni 



,-1 dq c 
dn f 



dn f 

The integrals are given by 



{h-h-h~e f )q- c 1 ^- (38) 



+ Wlhq-} dqCf 



k<kp 



and 



? \/-2 



(39) 



(40) 



for / = 2-4. From de/dd = 0, we obtain 



t/ c/ = - 2V 2 z / 2 4- 



-{h-h-he f )q c — 



(41) 



Here, note that while we take a(k) and ii as independent vari- 
ables for de/da(k) = and de/dd = 0, we take ny and a! as 
independent variables for the derivatives in eqs. (38) and (41). 
Equation (23) is rewritten using eqs. (37) and (40) as 



ny = - +h. 



(42) 



We solve eqs. (38), (41), and (42), and determine ey, d, and 
«/. By using eqs. (37), (39), and (40), we can rewrite eq. (36) 
as 



= h + n f e f + g - rifj l f -h- \V\h + Ucfd. 



(43) 



From the above equations, we find that the band structure of 
the conduction band is included only through the density of 
states in the present model, and a physical quantity, such as 
n c {k), depends on the momentum k only through e^. 

We can evaluate expectation values of physical quantities in 
the optimized wave function. For example, we obtain the mo- 
mentum distribution function n c (k) and n f(k) using eqs. (28) 
and (30), respectively. An important quantity is the jump 
An(kp) = An c (kp) + An f(kp) at the Fermi level; its inverse cor- 
responds to the mass enhancement factor. In the following, we 
call l/An(kp) the mass enhancement factor. 

In the next section, we show the calculated results for the 
model with a constant density of states. Before showing them, 
we discuss three characteristic regimes of the model, schemat- 
ically presented in Fig. 1, which do not depend on the details 
of the band structure. First, we consider a case with ey » W 
[Fig. 1(a)], where W is a typical energy scale of the conduc- 
tion band or half of the bandwidth in the next section. In this 
case, rif and the energy e per site is almost the same as 
the kinetic energy e^ n (n c ) per site of the free conduction band 
with n c = n. Second, we consider a case with ey <K —W, -U c f 
[Fig. 1(b)]. In this case, «/ - 1 and n c ^ n - 1. The energy is 



/ 



(b) Kondo regime 



E f » W 

n f ~0,e~e kin (n=n) 



/"•--•--•- 

e f «-W,-U cf 
rif~ 1, 



e ~ e kin( n c= n-\)+H f +(n-\)U cf 



(c) intermediate-valence regime 
c 

-Utf« -W 
n f ~ 2—n, e ~ (2—n)£ f 



Fig. 1. (Color online) Typical electron configurations in three character- 
istic regimes: (a) rif =s regime, (b) Kondo regime, and (c) intermediate- 
valence regime. ekin(«r) denotes the kinetic energy per site for the free con- 
duction band with n c . 



approximately given by e eidn(«c = n — 1) + e/ + (n - l)U c f- 
We call this regime the Kondo regime. For n/ — > 1, we obtain 
q — > 0, q c — > 1, q/ — > 0, and a(k-p) diverges as a(kp) ~ q~ l . 
By using them, we find that, for rif — > 1, An c (kp) — > 
and An j {kp) — > 0, that is, the mass enhancement factor be- 
comes large. This mass enhancement for rif — > lis consis- 
tent with the previous result for the PAM. Third, we consider 
a case with a moderate ey and a large t/ c y, more explicitly, 
-U c f « e/ « -W [Fig. 1(c)], In this case, / and conduc- 
tion electrons tend to avoid each other; thus, rif + n c /2 1 
and d - 0. That is, ny ^ 2 — n and n c ^ 2n — 2. Here, we 
call this regime the intermediate-valence regime. In this case, 
both / and conduction electrons are almost localized, and the 
energy is e (2 - n)ey. For ny + n c /2 — > 1 and d — * 0, 
we obtain q — > 0, q c — > 0, and gy — ♦ 0. By using them, we 
can show that the mass enhancement factor becomes large in 
this intermediate-valence regime. This mass enhancement in 
the intermediate-valence regime is not realized in the ordinary 
PAM and is a result of the effect of U c f. 

3. Results 

Now, we show our calculated results. Here, we consider a 
simple model for the kinetic energy: the density of states per 
spin is given by p(e) = 1/(2 W) for - W < e < W; otherwise, 
Pie) = 0. 

Figure 2(a) shows ny as a function of ey for several t/ c y 
values for V/W = 0.1 and n = 1.75. For a large U c f, we rec- 
ognize the three regimes mentioned above. A first-order phase 
transition occurs from the Kondo regime to the intermediate- 
valence regime or to the ny regime for U c f/W > 0.89. We 
observe hysteresis by increasing and decreasing ey across the 
first-order phase transition point, and here we show the values 
of the state that has the lower energy. Figure 2(b) shows the 
number of interacting electron pairs d through U c f per site. 
For a large I/ C y, the conduction and / electrons tend to avoid 
each other and d is suppressed. Figure 2(c) shows the valence 
susceptibility xv = —drtf/de/ as a function of ey. The va- 
lence susceptibility enhances around the boundaries of three 
regimes for a large U c f. For U c f = 0, such a boundary is 
not so clear. Figure 2(d) shows the mass enhancement factor 
l/An(kp) as a function of ey. In addition to the enhancement 
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Fig. 2. (Color online) et dependences of (a) nf, (b) d, (c) xw, and (d) 
1/An(foj) for V/W = 0.1 and n = 1.75. U cf /W = (dotted lines), 1 (dash- 
dotted lines), 2 (dashed lines), and 3 (solid lines). 



20 
15 



sfio 

< 



«=1.75 





U cf /W=0 
U cf /W=l -■- 
U cf /W=2 
U cf /W=3 



(l-« / /2)/(l-n / ) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

"/ 

Fig. 3. (Color online) 1/Arc(fcp) as a function of nt for V/W = 0.1 and 
n = 1.75. Ucf/W = (dotted lines), 1 (dash-dotted lines), 2 (dashed lines), 
and 3 (solid lines). The thin line is (1 - n//2)/(l - n/). The vertical line 
indicates re/ = 2 - n. 




(b) U cf IW=l 



(d) U cf /W=3 



-1 -0.5 0.5 1 -1 -0.5 0.5 1 



e k /W 



e k IW 



Fig. 4. (Color online) Momentum distribution functions n c {k) (dashed 
lines) and nAk) (solid lines) as functions of for V/W = 0.1, n = 1.75, 
and n f = 2 - n = 0.25. (a) U cf /W = 0, (b) U c f/W = 1, (c) U cf /W = 2, and 
(d) Ucf/W = 3. 




n=1.75 



0.5 1 1.5 



U cf /W 



for Hf — > 1 as in the ordinary PAM, we find another region, 
that is, the intermediate-valence regime nj 2 - «, in which 
the mass enhancement factor becomes large. This enhance- 
ment, in particular, a peak as a function of e/, is not expected 
for the PAM without U c f. Our theory may be relevant to the 
large effective mass in the intermediate-valence compounds 
<?-YbAlB4 and /?-YbAlB4 and the nonmonotonic variation in 
the effective mass underpressure in CeCu2Si2- 

To clearly observe the effect of U c j on the mass enhance- 
ment, we show l/An(£ F ) as a function of nf in Fig. 3. The 
thin line, which almost overlaps with the U c f = data, rep- 
resents the mass enhancement factor, given by eq. (1), i.e., 
(1 - n//2)/(l - n f ) obtained for the PAM with U cf = and 
g — 1 . Note that, in the present theory, g ± 1 even for U c f = 0. 
By increasing U c f, 1/An(kp) becomes large, particularly in 
the intermediate-valence regime «/ ^ 2 — n. 

In Fig. 4, we show the momentum distribution functions 
n c (k) and n f(k) for n — 1 .75 and n/ = 2- n = 0.25 for several 
values of U c f. For U c f - 0, the jump at the Fermi energy e^ F = 
(n - l)W = 0.75 W is much larger for rtf(kp) than for n c (kp), 
that is, the quasiparticle weight is mainly composed of the 



Fig. 5. (Color online) Inverse of peak value Xy (dashed line) of valence 
susceptibility [see Fig. 2(c)] and jump 5tif (solid line) in «/ [see Fig. 2(a)] at 
first-order phase transition as functions of U c f for V/W = 0.1 and n = 1.75. 
The circle represents the critical point. 

/-electron contribution. For a large U c f, the jump becomes 
small for both «/(&f) and n c (kp), and the mass enhancement 
factor becomes large, as shown in Fig. 3. 

Figure 5 shows how we determine the critical point of the 
valence transition. In this figure, we draw the inverse of the 
peak / ty eak of the valence susceptibility and the jump 6nj in 
nf at the first-order valence transition. Both of them should 
become zero at the critical point, and indeed, we find that they 
become zero at the same U c f. 

In Fig. 6(a), we show the valence susceptibility x\ as a 
function of £f and U c f for n = 1.75. In this figure, we also 
draw the first-order valence transition line and its critical 
point. The crossover lines, represented by the dotted lines, are 
determined by comparing the energies of the three extreme 
states: n/ = 0, «/ = 1, and n/ + n c /2 = 1 with d — 0. The 
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Fig. 6. (Color online) (a) XV and (b) 1/A«(&p) as functions of et and U c t 
for n = 1.75 with V/W = 0.1. The solid lines represent the first-order valence 
transition line. The solid circles denote the critical point of the valence tran- 
sition. The dotted lines indicate crossover lines determined by comparing the 
energies of the three extreme states (see text). 



crossover lines are given by 

e f = -(« - l)U cf + e kin (n c = ri)- ekinOc 
between the Kondo and rif regimes, by 

ekin(« f = «) 



2-n 



1), (44) 



(45) 



between the intermediate-valence and rif ^ regimes, and by 

e\dn(n c =71-1) 



e f = -Ucf - 



1 



(46) 



between the Kondo and intermediate-valence regimes. The 
crossover line between the intermediate-valence and rif — 
regimes does not depend on U c f, The other crossover lines 
are straight lines with finite slopes. Between the Kondo and 
rif =i regimes, the slope is -l/(n - 1) and does not depend 
on the band structure. Between the Kondo and intermediate- 
valence regimes, the slope is -1 independent of both the band 
structure and filling n? ] The region where xv becomes large 
is captured well by the crossover lines obtained by such a 
simple consideration. The first-order valence transition occurs 
only from the Kondo to intermediate-valence or to tit 
regimes within the U c f range presented here. Note that the 
valence transition can occur also between the intermediate- 
valence regime and the rif - regime for a smaller n. 20) 
Figure 6(b) shows the mass enhancement factor l/An(kp) as 
a function of ef and U c f. A large mass enhancement occurs 
in the intermediate-valence regime in addition to the Kondo 
regime. Here, note that the large mass enhancement occurs 
in the middle of the intermediate-valence regime. Thus, this 
enhancement is not due to valence fluctuations. In CeCu2Si2, 



1.20 



1.15 



1.10 



1.05 



1.00 




-1.5 
e f /W 



-0.5 



Fig. 7. (Color online) Ratio Snf/6d (solid line) of jumps [see Figs. 2(a) and 
2(b)] at the first-order phase transition and slope -5U c f/Sef (small circles) 
of first-order phase transition line (see Fig. 6) for V/W = 0.1 and n = 1.75. 
The large circle indicates the critical point. 



the effective mass has a peak before the superconducting tran- 
sition temperature becomes maximum under pressure, and 
which is consistent with our theory provided that the pairing 
interaction of superconductivity is mediated by the valence 
fluctuations. The situation will also be similar for CeCu2Ge2 
if we can subtract the contributions of magnetic fluctuations. 

Finally, to verify the consistency of the present theory, we 
check the Claudius-Clapeyron relation for the first-order va- 
lence transition. 7 ' 22) This relation is given by 

6n f 5U c f 
5d 6e 



(47) 



where 5d denotes the jump in d at the valence transition, and 
dU c f/def is the slope of the valence transition line. In Fig. 7, 
we show the values of the quantities on the left and right sides 
of eq. (47). We can clearly see that the Clausius-Clapeyron 
relation holds in the present theory. Note that the Clausius- 
Clapeyron relation also holds for the crossover lines men- 
tioned above. For example, rif 1 and d - n - 1 for the 
Kondo regime and rif ^ and d - for the rtf ^ regime, 
and then, between these two regimes, 6rif/5d l/(n- 1). It is 
the slope -5U c f/d£f for that crossover line given by eq. (44). 

4. Summary and Discussion 

We have studied the extended periodic Anderson model 
with U c f by Gutzwiller approximation. We have found that 
the three regimes, that is, the Kondo, intermediate-valence, 
and rif ^ regimes, are clearly defined for a large 
U c f. Then, we have found that, in the intermediate-valence 
regime, the effective mass is enhanced substantially. Accord- 
ing to the present theory, the large mass enhancement in the 
intermediate-valence regime indicates a large U c f- Thus, our 
theory provides helpful information for searching a supercon- 
ductor with valence-fluctuation-mediated pairing. 

In this study, we have not considered the possible insta- 
bility toward a spin-density-wave state and a charge-density- 
wave state. Such a state would be realized in a portion of the 
parameter space, particularly, in a lattice without geometric 
frustration. 819 ' The extension of the present theory to such 
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states is a future problem. In the present theory for a uniform 
state, the effect of a lattice structure is included only through 
the density of states of the conduction band. Thus, our results 
may change little for a frustrated lattice with a similar density 
of states even if we consider the possibility of the density- 
wave states. 

In our theory, we expand both the conduction- and f- 
electron states in the basis states in real space; thus, it will be 
possible to include the onsite correlation between conduction 
electrons and other short-range correlations. These extensions 
are future problems. 
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Appendix A: Approximation for Determinants 

In this appendix, we introduce approximations for determi- 
nants to evaluate expectation values in the variational wave 
function. Although most of them have been derived in ref. 10, 
we repeat them for the readers' convenience. 

We consider the state 



\k 1 ---k N ) = cl---cl\0), 



(A-l) 



where c' k denotes the creation operator of a spinless fermion 
with the momentum k t . From eq. (8), we obtain 

{k\ ■ ■ ■ k' N \kx ■ ■ ■ k N ) 

= J] det[<£ (r)] detfo>*(r)] (A . 2) 
M 

Then, we approximate each product of determinants by the 
average, that is, 



|det[^(r)]| 2 - -L, 

L<-JV 



and 



(A-3) 



(A4) 



det[^,(r)] detfo(r)] * 0, 
for {k} + {it'}. 

For the kinetic energy, we need to evaluate another type of 
determinant. We consider 



c r \k\ ■■■k N ) = ^ J det[<p k (r)]cr'\ri ■■■r N ) 

M 

= J] det (r ' ) [<p k (r)]\n---r N . l ), 



(A-5) 



where 



det (r,) [^(r)] 



VkM') 
tpkM') 



fkMO 



<Pk N ( r N-l) 



(A-6) 



Then, for r' + r" , 

(ki ■ ■ ■ k N \cl,Cr»\ki ■ ■ ■ k N ) 

= 2 det (r,) [^(r)]det (r " ) [^(r)]. 



(A-7) 



{r}j>r',r" 



On the other hand, by using the expansion 

Cr>=Y J <Pk'(r')c k ', (A-8) 

we obtain 

(k\ ■ ■ ■ k N \c^cAk\ ■■■k N ) 

= E fk'( r ">fk'(r"){ki ■ ■ ■ k N \c\,c k '\ki ■ ■ ■ k N ) 

k' " (A-9) 

N 

= ^^(r')^,(r"). 

i=i 

Then, we approximate the products of determinants in 
eq. (A-7) by their average: 

det (r,) [^(r)]det (r ">[^(r)] 

1 * (A- 10) 

Y.viy^kAr"). 



-iCn 



i=\ 



Appendix B: Evaluation of 2 II a 2 (k) with Restriction 

In the canonical ensemble for an free-electron system 
with dispersion s k at temperature 1//3, the hole distribution 
function is given by 

{c k c\)= Yj e'^^'/ZiN), (B-l) 
{if}** 

where Z(N) is the partition function. It should be equivalent 
to that in the grand canonical ensemble, 

1 + e P(e k -VL) ' ^ B ' 2 ^ 

where fi is the chemical potential, and thus, 

WW 

By putting Ek = - In a 2 (k),fi = 1, and ef 1 = q, we obtain 

-l 
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